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STATEWIDE LANDSAT INVENTORY OF CALIFORNIA FORESTS 
Abstract 

As mandated by the California Forestry Resources Assessment and 
Policy Act of 1977 (Assembly Bill 452) , the California Department of 
Forestry (CDF) was required to undertake an Initial forest land assess- 
ment for completion by July 1, 1979. A more detailed assessment and 
program plan lii to be completed by 1987. As part of the Initial assess- 
ment, Imagery from NASA's Landsat resource mapping satellite was util- 
ized. A cooperative mapping effort was agreed upon between the CDF, 
NASA's Ames Research Center (ARC), and the Jet Propulsion Laboratory 
(JPL), Pasadena, California, In August of 1978. The statewide mapping 
for the Initial assessment was completed by July 1, 1979 and represented 
considerable time and cost savings over traditional statewide mapping 
methods. Six forest cover categories were mapped, along with 10 general 
land cover classes. To map the state's 100 million acres, 1.6 acre 
mapping units were utilized. Map products were created. Standing 
forest acreage for the state was computed to be 26.8 million acres. 


InttoducUlon 

The Caltfprnia Forest Resources Assessment and Policy Act was 
created by the enactment of California Law AB452 in 1977 (State of 
California 1977) . This Act required the California Department of 
Forestry to develop an initial statewide forest assessment by July 1, 
3979, and a complete assessment by 1987. The assessment included an 
evaluation of existing forest data issues involved in the management of 
various forest resources, current policy, and the need for developing 
data collection and analysis systems. The initial report, "California's 
Forest Resources," was published in June, 1980 (California Department 
of Forestry 1980). Serious data limitations were found In resource 
issue after resource issue. One deficiency is the lack of a usable 
comprehensive data base. 

In California there are numerous detailed forest maps of limited 
areal coverage, Including ones by the U.S. Forest Service, and various 
timber company lands usually without public access. At least three 
maps have statewide coverage: the Forest Resources map (Weislander and 

Jensen 1945), Kuchler's Natural Vegetation of California map (in Barbour 

1977) , and the CALVEG map developed by the U.S. Forest Service (USFS 

1978) . Only the CALVEG map was used extensively in the initial assess- 
ment. The Natural Vegetation map is published at a scale of 1:1 million 
and displays "potential" natural vegetation. The CALVEG map was pro- 
duced through visual photo-interpretation of Landsat false color com- 
posite photograpblc prints. This was completed without precise geometric 
control of the mapping surface using an 800 acre minimum mapping unit. 
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A*\thougli CALVEG iB useful for developing policy Infomwcion, Its 
polygonal nature, large mapping unit£<, imprecise geometry, and at 
times questionable interpretation, make it questionable as a data base 
for future assessment of an operational analysis system. 

Use of Landsat satellite digital data in the construction of a 
suitable map base was then considered. Several advantages are found. 

The data is digital, thereby allowing computer processing, and providing 
discrete information for much smaller mapping units (down to 1.6 acres). 
It is possible to obtain complete state coverage in a single nine day 
period which can be updated frequently. In addition, previous studies 
have shown Landsat can provide useful resource information in a timely 
fashion. 

These factors resulted in CDF's decision to undertake statew/ide 
forest mapping using Landsat. This effort provided CDF an initial test 
of the suitability of this technology in fulfilling their resource data 
needs. A cooperative agreement to accomplish this task was formalized 
in June of 1978, between the California Department of Forestry, NASA's 
Ames Research Center, and the Jet Propulsion Laboratory (California 
Forest Resources Mapping Program Project Plan, 1978, available from 
NASA Ames Research Center, Technology Applications Branch, Moffett Field, 
California 94035), This project was to be completed by July 1, 1979. 
Final products vjere to be Landsat based digital and photographic land 
cover data for California. In addition, digital terrain data were to 
be derived from Defense Mapping Agency/National Cartographic Information 
Center digital elevation data. 
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Introduction to Landsat 


The first Landsat satellite, ERTS-1, was launcheu on July 23, 1972. 
Since then there have been two replacements, Landsat 2 (launched Jan- 
uary 23, 1975) and Landsat 3 (March 3, 1978). Today both Landsats 2 
and 3 are operational. 

The special-purpose Landsat satellites have been designed to mea- 
sure Earth resource features. Landsat imagery has been used in such 
fields as geology, urban planning, hydrology, mining, agriculture, and 
forestry (NCSL 1978) . 

Landsat is actually a satellite system and is normally referred to 
in the singular, even though more than one satellite may be operating 
at a time. 

Landsat operates at an altitude of 570 n* ml. Its orbit is 
designed so that repetitive ground coverage is obtained every -.ine days 
when more than one satellite is operating, as their orbits are comple- 
mentary. The satellite is in a "sun synchronous" orbit. At any loca- 
tion, images obtained on different days will be at the same local time 
of day (the middle to late morning in the northern temperate latitudes) 
(uses 1979). 

The heart of Landsat is the Multi-Spectral Scanner (MSS) . It 
obtains images roughly 100 n. mi. on a side. Each image is made of a 
grid of discrete picture resolution elements called "pixels." Each 
pixel represents a ground rectangle of roughly 1.1 acres. For each 
pixel, four spectral values are measured in four broad bands: one in 

the green (0.5 to 0.6 microns), one in the red (0.6 to 0.7 microns), 


k 



and two In the near Infrared (0.7 to 0.8 and 0.8 to 1.0 microns). 

These values describe the integrated solar reflectance intensities for 
all objects within the pixel (tJSGS 1979). 

The 100 n. mi. x 100 n» mi. images are known as "frames” or 
"scenes," not as photographs, as they are actually obtained in a digital 
electronic format that is computer compatible. As such, they are 
recorded In original form on magnetic tape and can be visually displayed 
on photographic film or color video monitors. This is done by display- 
ing the green image band as blue, the red band as green, and one of the 
infrared bands as red. The "false color" image has the same visual 
appearance as a color infrared photograph. 

There are. two ways to analyze Landsat imagery. It can be photo- 
interpreted in much the same way as one interprets a color infrared 
photograph, or it can undergo computer analysis (NCSL 1978) . In the 
first approach the image is analyzed using spectral and spatial pattern 
recognition. This is the approach used for the CAL7EG map. The second 
approach relies upon automated spectral analysis. Advantages and dis- 
advantages exist with both methods. The computer approach is presently 
handicapped by a lack of spatial pattern recognition capability. On 
the other hand, manual photo-interpretation is subjective, and inter- 
preted results are often inconsistent between analysts. Computer pro- 
cessing yields quantified results. Computers are also very fast and do 
not tire of processing voluminous quantities of data. This allows auto- 
mated analysis of each of the 7-1/2 million pixels contained in a Landsat 
image. 
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Each method benefits from the use of collateral Information, such 
as aerial photography and ground knowledge of the image area. In fact, 
Landsat data works best in harmony with other data structures, especially 
in cUgltal form. 

Landsat imagery can yield cost savings when the area to be analyzed 
is large, and the level of information required is broad and moderate. 
Large area ground surveys and aerial photographic mapping techniques 
are relatively expensive both for data acquisition and analysis. Landsat 
image analysis of large areas can be relatively Inexpensive as fewer data 
sets are required for area-wide coverage and spectral analysis is reason- 
ably efficient (NCSL 1978) , 

Landsat data are available from both the U.S. Geological Survey's 
Earth Resources Observation System Data Center, Sioux Falls, South Dakota, 
and the Canadian Centre for Remote Sensing, Ottawa, Ontario. 

Project Method 

The Landsat analysis for this project was composed of several com- 
ponents. First, the image data had to be selected. Then a digital 
mosaic of all the Landsat Images covering the state was produced (Fig. 1) . 
The data were then broken into image blocks for spectral analysis 
(Fig. 2). The final steps involved production of final output products 
and project documentation. 

One of the Initial steps was selection of appropriate Landsat 
Imagery for analysis. A digital mosaic of August 1976 scenes of Southern 
California had already been prepared at the Jet Propulsion Laboratory for 
the Bureau of Land Management (3LM) . This month ' s Imagery was good 
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BtatewidCf except for a small cloud-coverod area along the northern 
coast. Extending the August 197h southern California mosaic northward 
was chosen ns an economical method for producing suitable project data. 

The year 1976 was a drought year* and although tree mortality was not 
visible in August, overgrazing and very low reservoir levels were quite 
evident, An April, 1977, scene was used in place of the '^'loud-covered 
north coastal area. In total, 31 handsat scenes were used to produce 
tlie final mosaic, including some scenes from the BLM work. 

Because of the scene geometry, spaeecraft~induced anomalies, and 
the overlapping nature of adjoining handsat scenes, a continuous Landsat 
digital mosaic must be produced. This is done by accounting for geo- 
metric and radiometric differences between scenes and removing redundant 
data in overlapping areas. Projecting onto a map base of fixed geometry 
gives a continuous statewide image that can be processed as a single 
image rather than 31 parts. Various edge matching techniques are used 
to digitally blend images. In addition, through resampling image ele- 
ments, pixel size was changed from the original 57 m x 79 m to 80 m x 80 m. 
This larger, 1.6-acre, pixel size reduces the pixels needed to cover the 
state by 38%, thereby reducing computer processing time. 

The data still to be mosaicked is large (62 million pixels) . The 
complete mosaic could not be produced in a single block, so it was com- 
pleted in 5 pieces, each yielding a complete geographic sub-block (Fig. 3). 
The NASA Jet Propulsion Laboratory produced the mosaic (Development of the 
California Mosaic, 1978, Jet Propulsion Laboratory, Pasadena, California 
91103, unpublished document) . 
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Each block vas sent to Ames upon completion, starting from the 
already processed southern end of the state and moving norths The 
data were subsectioned into 1“ long by 1® let quadrangles to simplify 
downstream processing and applications work* This convenient size is 
equal to one-half of a standard U.S, Geological Survey It 250, 000 scale 
topographic series map. A Lambert conformal conic projection was used 
for the data. 

In addition to Landsat image data, digital elevation, slope, and 
aspect data wore provided for each quadrangle. JPL developed this 
information from Defense ^b^pplng Agency (DMA) digital terrain tapes. 

The terrain data were also registered to the map base. These data were 
useful during analysis as a means of differentiating vegetation types, 
as vegetation variations often relate directly to terrain variations 
(Strahler, Estes, and Maynard 1980). 

JPL also provided a translation equation relating latitude and 
longitude to pixel locations. This made it possible to r-zadlly relate 
map features to their corresponding pixels in the image and terrain data. 

Ames Research Center analyzed the image of each quadrangle. The 
"unsupervlsed classification" method was selected as the means of analy- 
sis. In this method, a clustering algorithm is used to develop spectral 
clusters. These clusters are defined by a spectral mean and variance in 
each of the 4 Landoat bands. The analyst controls the formation of 
clusters by defining the number desired and the minimum degree of separa- 
tion between each. The spectral values of all pixels are then examined 
in a four-dimensional space defined by the Landsat MSS bands. The algo- 
rithm will partition the spectral space so that the mean value of each 
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cluster centers on some congregation of values In the four-dimensional 
spectral space. Bach pixel can then be defined os belonging to a 
spectral group. Ideally feakurefiS fdcked out by each spectral cluster 
can be easily related to some desirable land cover type. Difficulties 
result, however, when a spectral class Includes 2 types, for example, 
picking out both residential and vineyard areas, This particular spec- 
tral class may pick out these featutes well, but data use is limited 
when thematic classes receive names such as resldentlal/vlneyard. This 
results in some specialized stratification to differentiate the 2 com- 
ponents, or simply labeling the spectral class either residential or 
vineyard, depending on which is the major component. 

This description is generalized. The actual method undertaken 
employed seypt-al variations (Fig, 4). First, the total area of each 
1" X 1® quadrangle was not clustered as a unit. In order to eliminate 
confusion between spectrally similar, but thematically opposed, land 
cover types, clustering was applied by ecological areas (Fig. 5) 

(Newland, Peterson, and Norman 1980). These areas were delineated by 
CDF staff based on Kuchler's Map of Potential Natural Vegetation (in 
Barbour 1977) . Thirty-one ecological zones (ecozones) were developed 
in this manner, As can be seen in Fig. 4, they extend across the Landsat 
1° X 1® quadrangle boundaries. Each ecozone covers an area of similar 
and broad plan communities. Clustering analysis was carried out on each 

* 

ecozone. 

As mentioned, several quadrangles are contained in each ecozone. 

The computer capability available did not support simultaneous cluster- 
ing of portions of multiple quadrangles. Some form of data assembly and 
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compaction was required. A process known as "weighting" was carried 
out (Newland, Peterson » and Norman 1980). As may already be evident i 
the development of spectral cluster mean and variance statistics can 
be obtained without reference to the spatial relationships between 
pixels. Only spectral values need be recorded. In weighting, each 
spectral value eni^ountered is recorded, and a counter keeps track of 
the rate of occurrence. In the process, inter-pixel spatial relation- 
ships are lost, but a great deal of data reduction results. This 
weighting process was used to interrogate the pixels within each quad- 
rangle containing a given scozone, extract the spectral values of the 
pixels in the ecozones, and record them in a computer file. When this 
was done, the 1“ x 1° image quadrangles could be tc-moved from the com- 
puter system, thus simplifying storage requirements. The weighted file 
developed would undergo clustering, normally for 35 to 45 spectral 
classes, this range having been decermln>?d optimum for the mapping task 
at hand. 

Each Landsat 1® x 1" quadrangle Siad the spectral class statistics 
reapplied to it once the clustering process was carried out on each 
ecozone. Each pix >1 within au ecozone was matched against the spectral 
statistics for that ecozone and assigned to a spectral class using a 
maximum likelihood decision algorithm. The resultant "classified" 

Itnage obtained by the interrogation of each pixel in this manner has 
1 piece of information per pixel — the spectral class to which it was 
most similar, 

This raised the question of exactly what features were to be mapped. 
Prior to full-scale processing, several quadrangles were put through the 
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weighting, clustering, and classification processes described »bove. 
these were analyzed against the context of an anticipated classification 
scheme to derive a legend of types and a definition in categorical form 
that could be achieved. Happing of the rest of the state was then made 
to conform to the legend (Table 1). 

Spectral class identification for each ecozone was undertaken once 
the block of image quadrangles covering that ecozone was classified. 

Each spectral class was examined by CDF and NASA to determine which of 
the land cover classes it best represented. This was done by examining 
each class on a color display, locating where it was occurring on lf-2 
color infrared photography (scales varied between 1 j 30,00() and 1}130,000), 
and photo-interpreting the land cover at those sites. Examination of 
the land cover picked out by a spectral class at several sites was 
usually sufficient to allow the analysts to extrapolate their identifi- 
cation across the whole ecozone. The CDF foresters were Instrumental in 
accomplishing spectral class Identification due to their field experi- 
ence in the areas mapped. 

Fitting an appropriate land cover name from the mapping legend to 
a spectral class was not always straightforward, Occasionally, spectral 
confusion resulted when 2 cover types were represented by a spectral class. 
When possible, these confusions were resolved by stratifying an ecozone 
with the digital elevation data and assigning the appropriate label 
within each elevation band. This cype of stratification was used to a 
limited degree, recognizing that a great deal more specificity could be 
obtained with more time and effortj in general, classes were labeled on 
a "best fit" or majority basis. 


After all labeling of the classes in a classified quadrangle was 
complete, the spectral classes were grouped into the 16 basic land cover 
types. County boundaries were digitally encoded and overlaid on the 
data. Acreage tabulations of the land cover types by county, and 
a color-coded photographic map of the state's resources, were then 
created . 

Results 

As a result of the project, acreage ta'^ulatlons for all land coyer 
types mapped were produced (Table 2). Photographic color prints were 
made at 1:1,000,000 and 1:250,000 scales for each quad, the latter scale 
being useful as it is the same scale as the USCS topographic map series. 
All information is also stored on computer tapes. 

The total cost for the project has been estimated to be less than 
$300,000 '<^lgure Includes Industry equivalent salaries), or less than 
0.3<; per acre. 

The project was established in June 1978, and completed In June 
19/9. The Image analysis at Ames Research Center by NASA and CDF staff 
was accomplished in the 8 months between November 1978 and May 1979. 

The output Is presently undergoing further analysis in areas of 
special interest by CDF and other state and local agencies. The map 
output will be used with other data from the CDF to help complete the 
program plan required by 1987. 
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Table 1. Glassification scheme* 


1. 

Alkali flats 


2. 

Bare rock 


3* 

Barren 

vegetation. 

4. 

Water 


5. 

Other 

Ice, snow, or clouds. 

6. 

Grassland 

>5% vegetative cover; <10% tree canopy. 

7. 

Open shrub 

5-24% vegetative cover is shrub species; <10% 
tree canopy. 

8. 

Brush 

25% or greater of vegetation is shrub species; 
<10% tree canopy. 

9. 

Conifer 

>25% conifer canopy closure; hardwoods <20% of 
the tree species present. 

10. 

Conifer /hardwood 

>25% canopy closure; conifers comprise 50-80% 
of the stand; hardwoods 20-40% of the stand. 

1 • 

Hardwood 

>25% hardwood canopy closure; conifers <20% of 
the tree species present. 

12. 

Hardwood / conlf er 

>25% canopy closure; hardwoods 50-80% of the 
stand; conifers 20-49% of the scand. 

13. 

Conifer /woodland 

10-25% conifer canopy closure; hardwoods <50% 
of tree species present. 

14. 

Hardwood /woodland 

10-25% hardwood canopy closure; conifers <50% 
of tree species present. 

15. 

Agriculture 


16, 

Urban 
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Table 2. Total acres by class by county 
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Lassen 1,221 34,413 3,226 401,269 199,121 1,241,319 501,126 
Los Angeles 4,413 149,931 20,145 562,418 15,396 15,847 404,258 588,383 791,278 
Madera 38,160 6,699 250,924 106,556 195,493 2 78,224 
Marin 48,108 10,976 19,657 16,264 133,635 12,272 13,812 
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San Joaquin 5,795 507,811 61,134 298,634 30,386 
San Luis Obispo 67,259 99,094 138,182 1,083,842 616,317 
San Mateo 46,923 1,286 54,468 13,036 63,157 6,899 11,222 
Santa Barbara 35,390 40,723 31,068 31,241 701,489 616,423 
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Total acres 3.5 1.9 6.4 1.7 2.0 4.7 17.1 23.4 13.1 

(In millions) 
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San Joaquin 1,077 289 63 3,830 909,019 
San Luis Obispo 119,022 40,436 17,085 2,053 925 8,715 2,192,930 
San Mateo 48,775 28,015 12,130 21,789 15,314 35,200 358,214 
Santa Barbara 170,499 60,473 1,501 14,109 43,349 1,746,265 


Table 2 . Concluded 
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Total acres 4.9 3.4 3.5 2.8 3.2 9.0 0.7 0.5 101.5 

(In millions) 
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Figure 1.- Landsat scenes covering California. 
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Figure .1° quadrangles of California. 
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Figure 5*- Ecozones q£ California. 
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